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We investigate the electron transport properties of a model magnetic molecule formed by two 
magnetic centers whose exchange coupling can be altered with a longitudinal electric field. In 
general we find a negative differential conductance at low temperatures originating from the different 
scattering amplitudes of the singlet and triplet states. More interestingly, when the molecule is 
strongly coupled to the leads and the potential drop at the magnetic centers is only weakly dependent 
on the magnetic configuration, we find that there is a critical voltage Vc at which the current becomes 
independent of the temperature. This corresponds to a peak in the low temperature current noise. 
In such limit we demonstrate that the quadratic current fluctuations are proportional to the product 
between the conductance fluctuations and the temperature. 



An intriguing aspect of electronic transport is the inter- 
action between the current electrons and the internal de- 
grees of freedom of the conductor. Atomic positions and 
vibrations are certainly at research center-stage, electro- 
migration being the most obvious example of interplay 
between the current and the atoms motion. The sit- 
uation becomes even more intriguing at the nanoscale, 
where quantized vibrations can be detected by measur- 
ing the electron current and its derivatives with respect to 
the applied bias. This is the principle of inelastic electron 
tunneling spectroscopy (IETS). Furthermore also the re- 
verse effect is possible, namely one can control the atomic 
positions of a nano-object by exciting appropriately some 
vibrational modes. Current-induced chemical reactions^ 
and nano-catalysi:^ on surfaces are among the most ap- 
pealing potential applications of this field. 

Equally important is the interplay between the electron 
current and the magnetic texture of a magnetic device. 
Such an interplay underpins the giant magnetoresistance 
effectP and its reverse, i.e. current induced magnetiza- 
tion dynamics^. Considerably less investigated are the 
same phenomena at the atomic scale. This is mainly due 
to the intrinsic difficulties of both manipulating and de- 
tecting a few spins. In addition, magnetic excitations 
occur at energies lower than those involved in molecular 
vibrations, so that the measuring temperatures are often 
rather low. Still there are notable examples, such as ul- 
tra low temperature IETS of magnetic atoms on surfaces^ 
and of two-probe devices incorporating single magnetic 
molecules^. 

A new exciting prospect for scaling down spin- 
dynamics to the atomic level may be given by the ability 
of manipulating the magnetic configuration of a molecule 
with an electric potential instead of an electric cur- 
rent. Electrically induced alteration of the exchange cou- 
pling has been already predicted for two-centers magnetic 
molecules^ and nanowires^, and it is essentially based on 
the fact that the Stark shift of a magnetic object may de- 
pend on its magnetic state. This effect can be a crucial 
ingredient for the physical implementation of quantum 
computing based on spin el 9 * 10 ! 

An intriguing question is whether or not the depen- 
dance of the exchange coupling over an electrical poten- 



tial in a magnetic molecule can be detected electrically. 
This is the goal of our letter where we investigate the 
current- voltage, I-V, curve of a two-terminal device in- 
corporating a two-center magnetic molecule in which the 
exchange coupling changes with bias. Importantly we 
find that, in particular conditions of coupling between 
the molecule and the electrodes, there is a critical volt- 
age Vc at which the current becomes independent of the 
temperature. This is accompanied by a negative differen- 
tial conductance (NDC) at low temperature originating 
from the difference in scattering amplitude of the differ- 
ent spin-states of the molecule. 




FIG. 1: (Color on line) The model system investigated: a 
dimer of magnetic atoms (red), carrying respectively spin Si 
and 5*2, is attached to two ID non-magnetic electrodes (light 
blue). The Hamiltonian for the dimer is Hdim ■ The scattering 
region (dashed box) includes the dimer and six atoms of the 
electrodes and it is described by the Hamiltonian matrix Hs • 

In figure [T] we show the simple model system inves- 
tigated, which comprises a di-atomic magnetic molecule 
sandwiched between two one-dimensional non-magnetic 
electr odes. The system is described by the s-d 
modePilll where spin of the current carrying s-electrons 
is exchange-coupled to the local spins Si and S2 (d) of the 
two atoms in the dimer. Si and S2 are treated as classical 
variables and their orientation determines the scattering 
potential for the s-electrons. These are described by a 
tight-binding Hamiltonian with a single s-orbital per site 
at half-filling. The on-site energy and hopping integral 
in the electrodes are £l = 2 eV and 7l = — 2 eV, a choice 
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which maintains the system far from the van Hoof sin- 
gularities at any voltage investigated. No local spins are 
present in the electrodes so that their electronic structure 
is not spin-polarized. The Hamiltonian of the dimer is 
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where hf^ is the on-site Hamiltonian matrix of the i-th 
atom of the dimer, 7 is the hopping parameter and 
(cf ) is the creation (annihilation) operator for an elec- 
tron with spin a (T,4) at the site i. We have defined 
h T = I e o + U{pi- Po)]S a p - J sd Si ■ (<r) a p, where a are 
the Pauli matrices, pi is the total occupation of i-th site, 
po = 1 is the site occupation in the neutral configuration, 
U = 1 eV is the atomic charging energy and J s( j = 2 eV is 
the exchange parameter between the s-electrons and the 
local spins. In our calculations we consider eo = 2 eV, 
7 = —0.1 eV, = 1. In absence of spin-orbit interac- 
tion and spin-polarization of the electrodes the scattering 
potential is determined only by the mutual angle, 9, be- 
tween the two local spins. Finally, the dimer and the 
electrodes are coupled by the hopping integral 7c. In 
particular we explore the two cases in which 7c = 1 /47l 
and 7c = 1/2tl. These parameters are only illustrative 
and have been chosen in order to maximize the differ- 
ence in conductance between different spin-states of the 
molecule (6 — vs 9 = n) . 

The non-equilibrium Green function method^ applied 
to our tight-binding HamiltonianfSi is used to calculate 
the transport properties. The central quantity is the re- 
tarded Green's function of the scattering region 



G{E)= lim [{E + i v ) - #dim 
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where E is the energy, Hg is the Hamiltonian matrix 
of the scattering region and El (£r) is the self-energy 
of the left- (right-) hand side electrode. This latter 
describes the interaction between the scattering region, 
which includes the dimer and 6 atoms of the electrodes 
(see Fig. [IJ, and the electrodes. G(E) enters in a self- 
consistent procedure to evaluate the stationary occupa- 
tion of the scattering region and once convergence is 
achieved the two-probe microscopic current, i(V), at the 
voltage V is extracted from the Landauer formula!^. 

Since at any given temperature the angle between the 
magnetic moments in the dimer fluctuates, for any mi- 
croscopic quantity q we can define its macroscopic coun- 
terpart, Q, as the thermal average over all the possible 
angles 



Q(V) = (q(V)) 



d£ 



(3) 



thus that if q = i one obtains the macroscopic current, 
/. Here £ 12 is the dimer magnetic energy, which writes 



and £ m i n (£ max ) is its minimum (maximum) value. In the 
equations Q above Jdd is the exchange energy between 
the two spins, which in turns is a quadratic function of 
the electrical potential difference between them, Vd- This 
latter is an intrinsic function of both V and 9. Finally the 
constants a and b are fixed to the values of a — 0.001 eV 
and b = —0.8 eV/V 2 . Note that the functional depen- 
dance of Jdd over Vd implies a critical voltage at which 
the exchange energy changes sign, i.e. the magnetic cou- 
pling turns from ferromagnetic to antiferromagneticP^. 

We begin our analysis by investigating the microscopic 
quantities, i.e. the current i(V,9) and the dimer internal 
potential drop Vd- In Fig. [2] we present, for both choices 
of coupling 7c, i-6 for different voltages and Vd-V for 
different angles 9. In the case of 7c = 1/47l the cur- 
rent varies as i(V,9) ~ [ia + iiCOs9]V, with io and i\ 
two constants. At the same time Vd is only weakly de- 
pendent on the internal spin configuration, i.e. the Vd-V 
curve changes little with the angle 9 [Fig. [2jb)] . In con- 
trast for 7c = 1/2tl the current peaks at approximately 
6 = 7r/2 with both the parallel and antiparallcl config- 
urations being low conducting. Again the amplitude of 
the current variation over 9 increases with bias, although 
only moderately in this case. Furthermore for this situ- 
ation Vd, which is still linear with V, is rather sensitive 
to the angle between the two spins. These differences 
affect dramatically the macroscopic current, /, that we 
calculate next. 
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FIG. 2: Microscopic transport quantities. The microscopic 
current, i, as a function of the angle, 8, is shown in panels 
(a) and (c) for different voltages V. Panel (b) and (d) show 
the internal potential drop, Vd, as a function of the external 
bias and for different angles. Panels (a) and (b) are for 7c = 
1/4tl, while (c) and (d) are for 7c = 1/2tl. 



The macroscopic I-V curves for the two cases are pre- 
sented in the panels (a) and (c) of Fig. [3] while the 
panels (b) and (d) report the current quadratic fluctu- 
ations AI — (i 2 ) — I 2 still as a function of bias. In all 
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cases we study the electrical response in the temperature 
range 1-15 K. The most interesting behaviour is found 
for 7c = 1/4tl, from which we start our discussion. 





0.07 




0.06 




0.05 




0.04 


I 


0.03 






0.02 




0.01 









0.25 




0.2 




0.15 








0.1 








0.05 




°, 



- AW 

: M 


1 1 ____ 


/*•?/ 

- P 


V z 


7 , i , i 


c 

1 



1 1 1 






1 K 


/ 


3 K 

5 K 




7K 


S 1 1 


15 K 



0.1 0.2 0.3 

V(Volt) 




0.1 0.2 0.3 

V(Volt) 



FIG. 3: Macroscopic transport quantities. The I-V curves 
are presented in panels (a) and (c), while panels (b) and (d) 
display the current quadratic fluctuations AI = \J (i 2 ) — I 2 
as a function of bias V . Panels (a) and (b) are for 7c = 
1/4tl, while (c) and (d) are for 7c = I/27L. Note that in (a) 
there is a critical voltage Vc at which the current becomes 
independent of the temperature. 

Figure |3ja) reveals two remarkable features. First we 
note that there is a pronounced NDC at about 0.2 Volt, 
which is well evident at 1 K, it weakens as the temper- 
ature increases and finally disappears at 15 K. Interest- 
ingly the scaling of the electrical current with the temper- 
ature is opposite at the two sides of the NDC: it decreases 
as the temperature is enhanced before the NDC while it 
grows with T for voltages just after the NDC. The same 
NDC is present also in the case of stronger coupling with 
the leads [7c = 1/27l, Fig- [3jc)] at the somewhat lower 
voltage of about 0.1 Volt. In this case however the NDC 
is much less pronounced and disappears already at 3 K. 

The second and most striking feature of Fig. [3ja) is 
the presence of a critical voltage, Vc, at which the current 
becomes independent of the temperature. Such a voltage 
is in the vicinity of the NDC and correlates well with the 
peak in the current quadratic fluctuations [Fig. [3jh)] at 
low temperature. Note that this second feature is absent 
in the case of strong coupling to the leads. 

All these aspects can be easily understood by relating 
the microscopic quantities of Fig. [2] with the average of 
equation pi). Let us consider the case of 7c = 1/4tl 
first. In general the macroscopic current I(V) is deter- 
mined by the microscopic currents i(V, 9) of those con- 
figurations in which the system spends most of the time. 
The equations Q tell us that the ferromagnetic configu- 
ration is energetically favorable at low bias, while it is the 
antiferromagnetic to dominate at higher voltages (for 



larger than ±y/—a/b). This means that as the external 
bias increases the average current becomes progressively 
dominated by antiferromagnetic configurations to the ex- 
penses of the ferromagnetic ones. Since the microscopic 
current for 9 = 7r is always considerably smaller than 
that for 9 = [see Fig. |2ja)], this results in a decrease 
of the macroscopic current as a function of bias, i.e. in 
the NDC. Note that this particular NDC is not of mi- 
croscopic electronic origin since the microscopic currents 
i(V, 9) are monotonic in V for every 9. 

The fact that the exchange coupling changes sign as 
a function of the bias produces the second important 
feature in the macroscopic I-V curve. In fact when 
the potential drop between the two magnetic atoms is 
U(j = ±y/—a/b, then the parallel and antiparallel config- 
urations of the magnetic molecule become energetically 
degenerate. This means that now no magnetic energy 
scale enters into the problem and the system spends an 
equal amount of time in any spin configurations regard- 
less of the temperature. In general is proportional to 
the external bias V. Therfore one expects the existence 
of a universal external bias Vc = V[£i2(vd) = 0] such 
that £12 = and the macroscopic currents becomes in- 
dependent of the temperature, as indeed demonstrated 
in Fig. (3ja). However there is a second condition for this 
to happen, i.e. should be independent of the angle 9. 
This is not satisfied for 7c = 1/27l [see Fig. [2jd)] and 
as a consequence the I-V curves remain temperature de- 
pendent at any bias. 

From our discussion it is now clear that if is pro- 
portional to V and weakly dependent on 9, then there 
will be a critical voltage Vc at which any macroscopic 
quantity becomes temperature independent. Figure [3]( a) 
illustrates this feature for the current and the same is 
demonstrated in Fig. |4^b) for the conductance G = ( ) . 
Interestingly one can also adopt a different definition for 
the macroscopic conductance, namely that of the bias 
derivative of the macroscopic current dl/dV. Such a 
quantity is presented in Fig. |4^a) and as expected it ap- 
pears sensibly different from G. Interestingly both G and 
dl/dV are, in principle, accessible from experiments, and 
one may wonder whether some general conclusions can be 
taken by measuring the two quantities independently. 

In general, by taking the equation ^ and formally 
deriving Q with respect to the bias V we find 

dq 

W 



dQ 
dV 



/crT = 



d£ Vd 



dV 



-Q 



dV 
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where fee is the Boltzman constant. If one now considers 
q = i then the equation (|5| establishes a general relation 
between the conductance fluctuations and the correlation 
function between the current and the magnetic energy. 
Such a relation is drastically simplified when the mi- 
croscopic current has the typical spin-valve dependence 
i(V, 9) ~ [io + i\ co&9]V and «d is linear with V . In this 
situation (encountered here for 7c = 1/4tl) one finds 

or 
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FIG. 4: Differential conductance calculated as (a) the deriva- 
tive of the macroscopic current, I, with respect to the bias 
or as the (b) thermal average, G, of the microscopic conduc- 
tance. In panel (c) we show the product of (G — 81/ dV) 
with the temperature T, as a function of bias. Results are 
presented for 7c = 1/4tl. Note that [G — 8I/dV]T is pro- 
portional to the square of the current quadratic fluctuations 
Al 2 [panel (d)] of figure 3(b). 



i.e. that the conductance fluctuations rescaled by the 
temperature are proportional to the squared current fluc- 
tuations. A numerical proof of such a relation is provided 
in the panels (c) and (d) of figure [4] 



In conclusion we have investigated the temperature- 
dependent electronic transport through a model diatomic 
magnetic molecule, in which the exchange coupling be- 
tween the two magnetic centers is a function of the bias. 
This presents two remarkable characteristics. First, if the 
potential drop between the two magnetic centers is only 
weakly dependent on the angle between their magnetic 
moments and it is linear in V , then there is a critical volt- 
age Vc at which the macroscopic current becomes tem- 
perature independent. Secondly, if in addition the micro- 
scopic current has a form i(V, 9) ~ [io + i\ cos 6]V, then 
there is a universal relation between the temperature- 
rescaled conductance fluctuations and the quadratic cur- 
rent fluctuations. Both these effects are a unique finger- 
print of the dependance of the magnetic energy upon an 
external bias and can be used as a tool for detecting such 
a dependence. 
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